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Dynamical roughening of interfaces has received much attention in recent years. 
However, experiments have been restricted to one dimensional (Id) systems. Moreover, 
theoretical studies of the two dimensional (2d) case have been highly inconclusive. Here 
we introduce an experimental 2d system, with which the theories can be tested. As is 
shown, the surface of a 2d pile of rice shows roughening behaviour in both space and 
time, with a roughness exponent a.2d = 0.39(3) and a growth exponent fi^d = 0.27(3). 



1.1 Introduction 

Roughening phenomena of interfaces have been studied extensively in recent 
years due to their wide range of applicability. Rough interfaces appear in such 
diverse systems as flux propagation in superconductors [Q , the burning of papers 
Q , diffusion waves || , bacterial colonies Q| , flow through porous media Q and 
many more Q . Even though all of these systems have very different microscopic 
physics governing the processes, they can be described by simple models from 
a very small number of universality classes ||. The most famous such model is 
described by a non-linear diffusion equation known as the Kardar-Parisi-Zhang 
(KPZ) equation §: 



d t h(x, t) = vAh(x, t) + \(X7h(x, t)f + rj(x, t), (1.1) 
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where v is the diffusion coefficient, 77 is a noise term, A quantifies the non-linearity 
and h{x, t) is the position of the interface. In one dimension, the scaling behavior 
of an interface governed by the KPZ equation can be analytically solved. The 
roughness is parameterized by the width of the interface given by: 

a(t,L) = {((h(x,t)-{h(t)) L f) L ) 1/2 . (1.2) 

Here, (-)l denotes the average over the interface in space. For a self-affine 
surface, the width growth as a power law in time a ~ t 13 , until saturation is 
reached when the correlation length becomes comparable to the system size || . 
This growth exponent, (3, characterizes the dynamics of the process. After the 
saturation time, the width is constant in time at a value a sat (L) ~ L a , which 
grows as a power law with the system size 0. This roughness exponent, a, 
characterizes the structure of the interface. For the KPZ equation, one obtains 
a = 1/2 and f3 = 1/3 §. 

For a multi-dimensional KPZ system however, the theoretical situation is un- 
clear. Analytical treatments of the KPZ equation only exist in approximations 

^, Eo], and results from numerical simulations vary greatly ||. Similarly, 
experiments have up to now been restricted to a single dimension. The experi- 
mental problem of a two dimensional rough surface asks for a surface reconstruc- 
tion technique with enough spatial resolution to span some orders of magnitude, 
while at the same time having the temporal resolution to capture the dynamics 
of the process, which is not easily achieved. Secondly, a system has to be found 
that exhibits KPZ roughening in 2d and is accessible experimentally. Presently, 
there is some interest to combine the exact results on the KPZ equation with 
the concept of self-organized criticality (SOC) fTl|| . The surface of a sandpile, 
which is the archetypal system to study SOC, can be mapped onto a system 
which follows KPZ dynamics p2[ . This is intriguing since it brings together two 
established fields of research, however has not yet been tested experimentally. 
We study here the front of a 2d rice-pile and its roughening behavior, showing 
that it does indeed obey KPZ dynamics. We choose rice, since it has been shown 
in Id that a rice-pile does indeed show SOC behaviour jlq] . With the surface of 
a rice-pile established as a roughening system, we can extend the study further 
to include the full 2d surface of the pile spanning an area of ^lxlm 2 . Using this 
system, we can determine roughening and growth exponents in 2d and compare 
them with theoretical predictions ||, ||, [l0| . 

In section 2, we will discuss the experimental setup including the surface 
reconstruction technique based on active-light stereoscopy, as well as the grow- 
ing mechanism of the pile. In section 3, we develop the analysis techniques, 
with special emphasis on the generalization from known Id methods to the 2d 
problem. In section 4 we present the results of the rice-pile experiment. There, 
it is first shown that the front behaviour in Id does in fact obey KPZ scaling 
before discussing the 2d results. Those results are used to put constraints on 
theoretical results for 2d KPZ behavior. 
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1.2 Experimental setup 

The rice-pile is grown by dropping rice, uniformly distributed along a line using 
a custom-built dispenser. The dispenser consists of a distribution board and 
a sowing machine. In the sowing machine, an eccentric rotor keeps the rice in 
motion such that a steady flow of rice is achieved at the rate of ^5 g/s. This 
flow of rice is subsequently distribute d al ong a line of 1 m in the distribution 
board using simple geometry (see Fig. [O]). The principle is related to that of a 
pin-board producing a Gaussian distribution. 
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Figure 1.1: A schematic image of the distribution board. Rice is dropped from a 
single point on the top and subsequently divided into even compartments. At the end 
a line of rice uniformly distributed in 64 intervals is obtained, which is used to grow 
the rice-pile at a rate of ~5 g/s. 



In order to study the surface properties of the rice-pile, a 3d reconstruction 
technique was developed, based on active-light stereoscopy [Q. A set of colored 
lines is projected onto the pile at approximately right angles using an overhead 
projector. In the stereoscopic view the projector takes the place of the 
second camera passing its information to the camera via the colored lines. The 
camera itself is placed at an angle of 45 degrees to the surface of the pile and 
the projected lines. From this view-point the projected lines are deformed and 
can be used to determine the 3d structure of the surface in the same way as 
iso-hcight-lines do on a map. An example of such a reconstruction is shown in 



Fig. 1.2. Measurements on test objects show that a surface of lxl m can be 



reconstructed with an accuracy of 1-2 mm, which is comparable to the size of 
the rice grains and thus suited for the present purpose. The use of differently 
colored lines allows for better filtering and thus for better identification of the 
lines in the computer. 



1.3 Analysis methods 

As noted in the introduction, rough surfaces are often analyzed using the width 
of the interface to characterize its structure and dynamics. However, in order to 
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Figure 1.2: Reconstruction of the surface of a rice-pile. The white line indicates the 
position of the growing front. 



obtain reliable results, many experiments have to be averaged over, using such a 
method A more promising way of analysis, which has been extensively used 
in the analysis of Id experiments is via the two-point correlation function ra] 



C(x, t) = {(h(t t) - h& + x,t + t)) 2 ) 1 ^. 



(1.3) 



In both space and time the scaling behaviour of the correlation function is the 
same as that of the width thus making it possible to determine the growth 
and roughness exponents from C(x, t). In addition, the growth exponent can be 
determined from data obtained after the saturation time, since in the correlation 
function only time differences are important. 

When generalizing the method to 2d, computational difficulties arise. Be- 
cause of the number of points to compare, the number of operations to be carried 
out to determine the correlation function grows with the fourth power of the size 
of the surface. However, tests on small surfaces indicate that the radial aver- 
age of C(x,y,t) scales like the 2d local width, but due to the computational 
inefficiency we were using yet another method to determine the roughness and 
growth exponents. 

The power spectrum, or structure function | ]T6| , can be determined easily for 
Id and 2d systems from the square of the Fourier transform h(k x , k y ) of the local 
height h(x, y) 

S [k x i ky 

) = \h(k x ,k y )\ 2 . (1.4) 

Here, the computational load is just given by determining the Fourier transforms, 
which also in 2d only grows with the square of the size of the surface thus making 
it feasible to calculate the distribution function of the whole rice-pile surface. 
The square root of the integral of S (k x , k v ) over k-space, the distribution function 
cr(k x , ky), 

<r(k x ,ky) = I / S(k x ,k y )dk x dky) (1.5) 
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is equal to the rms- width of the interface [[t6| . Thus the distribution function also 
has the same scaling behaviour as the width and can therefore be used to deter- 
mine the roughness and growth exponents. Again, a radial average of a(k x , k y ), 
<r(fc), scales like the 2d local width, thus making comparisons with previous 
simulations of ballistic deposition models possible. Moreover, real 2d measures 
like the distribution function can also give information about anisotropies of the 
scaling in the x- and y-directions. 

The distribution function is also useful in investigations of the dynamics of 
the processes. In that case, the square of the Fourier transform h(u>) of the time 
dependence h(t) has to be determined. The Fourier transforms are determined 
using an FFT algorithm after padding the data with zeros to the next power of 
two. 



1.4 Results and discussion 

In order to determine that the rice-pile surface does in fact follow KPZ behaviour 
we first determined the roughness and growth exponents of the front of the pile, 
given by the line of equal height of the pile at 0.1 m. The distribution func- 




Figure 1.3: The behavior of the front of a propagating rice-pile. Both the correlation 
function (open symbols) and the distribution function (full symbols) show scaling in 
space and time over two decades. The resulting roughness and growth exponents are 
in excellent agreement with the KPZ universality class. 



tions determined in both space, cr(fc), and time, cr(uj), as well as the correlation 
functions C(x,t) can be seen in Fig. 1.3, where the values of a — 0.48(3) and 
(3 = 0.33(3) can be inferred. These values are in excellent agreement with the 
expectations from the KPZ equation thus establishing that KPZ behaviour does 
appear in SOC systems. 

The 2d distribution function, cr(k x ,k y ), which characterizes the roughening 
of the whole surface is shown in Fig. |1.4| on a triple- logarithmic plot. In the 
insert, the angular dependence of a power-law fit to a is shown. This indi- 
cates a dependence of the roughness exponent a on the direction, which shows 
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the anisotropy of the system. Such an anisotropy most probably arises from 
the growth mechanism of the pile, which is seeded from a horizontal line, thus 
breaking the symmetry of the x- and y-directions. It should be noted here that 
the exponents corresponding to the x- and y-directions do not have to agree with 
those determined in a Id analysis. This is because a(k x , 0) already includes data 
from the y-direction due to the complex nature of the Fourier transform. Thus, 
a{k x , 0) already presents an effective 2d measure. 




Figure 1.4: The 2d distribution function for the rice-pile surface on a triple- 
logarithmic plot. From a radial average, the roughness exponent can be determined. 
The insert shows an angular dependence of the roughness exponent, with an anisotropy 
in the x- and y-directions. 

The radial average of the distribution function, a(k), is shown in Fig. |l.5| a, 
where the value of the roughness exponent can also be determined. We ob- 
tain ct2d — 0.39(3), which is also in agreement with the average of the expo- 
nents determined as a function of angle. In addition, the temporal behavior, 



a{ui) is shown in Fig. 1.5b, where we determine the 2d growth exponent to be 
/?2d = 0.27(3). Both the roughness and growth exponents determined experimen- 
tally are in very good agreement with the conjecture derived from solid-on-solid 
models by Kim and Kosterlitz |l7| for higher dimensional exponents given by a 
= 2/(d+3) and (3 — l/(d+2). Numerical results from integrating the 2d KPZ 
equation vary greatly, with values of a^d ranging from 0.18 M to 0.39 (|] and 
/32d ranging from 0.1 jl(| to 0.25 Our experimental results are in good agree- 
ment with the numerical values of Amar and Family Q , as well as Bouchaud and 
Cates ||] corresponding to the high range of the values, while excluding most of 
the other numerical investigations into 2d KPZ behavior. 



1.5 Conclusions 

We have presented an experimental study on roughening in a 2d system. The 
surface of a rice-pile is measured with a reconstruction technique based on active- 
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Figure 1.5: (a) The radial average of the 2d distribution function, allowing the de- 
termination of the roughness exponent in 2d to be aid = 0.39(3) form a scaling regime 
spanning a decade and a half, (b) The distribution function in time, leading to a 
growth exponent of fod = 0.27(3). 



light stereoscopy. In Id, the fronts of the rice as the pile is grown shows excellent 
agreement with the Id KPZ universality class with exponents a = 0.48(3) and (3 
= 0.33(3) from a scaling-regime spanning more than two decades. Thus having 
established the KPZ nature of the system under study, we analyze the full 2d 
surface of the pile, where find a roughness exponent of (%2d = 0.39(3) and @2d 
— 0.27(3). This is consistent with numerical simulations for ballistic deposition 
models jl8[ [l^] and puts a strong experimental constraint on the available results 
on 2d KPZ simulations. Our results are in good agreement though with the 
results of Amar and Family ||, as well Bochaud and Cates || from numerical 
integration of the 2d KPZ equation. In addition however, we have studied the 
dependence of the exponent on the direction, where we find that the system is 
anisotropic with a somewhat higher exponent along the front direction. This is 
probably related to the difference between the two directions due to the growth 
mechanism. 
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